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ABSTRACT 

We present data from mid-infrared Keck Telescope imaging of 18 radio- 
selected ultra-compact H II region candidates at diffraction-limited resolution 1 . 
The goal of these observations is to determine the sizes, luminosities, and mor- 
phologies of the mid-infrared emitting dust surrounding the stellar sources. All 18 
sources were imaged at 11.7 fj,m and at 17.65 /im, and 10 of them were imaged also 
at 24.5 fim. All the sources were resolved. We have generated dust temperature 
and optical depth maps and combine them with radial velocity measurements 
and radio data (1.4 and 5 GHz) to constrain the properties of these star-forming 
regions. Half of our objects are excited by B-stars, and all our objects have 
derived types that are later than an 06 star. We find a significant correlation 
between infrared and radio flux densities, and a weaker one between infrared 
diameters and the central source ionizing photon rates. This latter correlation 



1 The data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, the University of California and the National 
Aeronautics and Space Administration. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 
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suggests that the more compact sources result from later spectral types rather 
than young age. Our new data may suggest a revision to infrared color selection 
criteria of ultra-compact H II regions at resolutions < 1". These 18 sources are 
part of a sample of 687 sources dominated by ultra-compact H II regions selected 
by matching radio and infrared maps of the first Galactic quadrant by Giveon 
and coworkers. The new mid-infrared images constitute a significant improve- 
ment in resolving sub-structure at these wavelengths. If applied to all of this 
sample our analysis will improve our understanding of embedded star-formation 
in the Galaxy. 

Subject headings: Galaxy: general — H II regions — infrared: ISM — ISM: dust 

1. Introduction 

Embedded star-forming regions constitute an important contribution to the overall 
Galactic star-formation (Wood & Churchwell 1989a; Churchwell 2002). The energetic ra- 
diation fields of massive stars ionize the surrounding interstellar medium, producing radio 
free-free radiation, and heat their natal dust cocoons, generating mid-infrared (MIR) radi- 
ation. Measuring the radio and IR luminosities of H II regions can therefore constrain the 
properties of the stellar objects embedded in them. The physical basis for the radio and 
IR emissions determines the spectral range most suited for these observations. The shape 
of the thermal free-free spectrum in the radio is such that the optical depth decreases at 
higher frequencies, enabling more sources to be detected at these frequencies. In the IR, the 
longer the wavelength, the cooler the objects that can be observed, making the observations 
more sensitive to more embedded sources in the mid and far-IR (FIR). Understanding the 
processes and the various evolutionary phases of embedded high-mass stars (e.g., protostars, 
young stellar objects, and ultra-compact [UC] H II regions) requires the study of large, com- 
plete samples. Statistically complete samples can establish the distribution of properties 
such as IR colors (e.g., Wood & Churchwell 1989b), radio spectral indices, spatial distribu- 
tions (e.g., Zoonematkermani et al. 1990; Becker et al. 1994; Giveon et al. 2005a; White, 
Becker, & Helfand 2005), and ultimately, the spectral types and the initial mass-function 
and formation-rate of massive stars (e.g., Wood & Churchwell 1989b). 

Numerous studies have addressed the issue of determining the stellar content, the evo- 
lutionary phase, and the dust properties of individual embedded sources or small to medium 
(< 100 sources) samples of them, by combining IR and radio observations: embedded sources 
and molecular clouds (Brand & Wouterloot 1991; Mookerjea et al. 1999; Saito et al. 2006), 
protostars and star-forming regions (Homeier 2003; Kraemer et al. 2003; Klein et al. 2005; 
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van der Tak, Tuthill, & Danchi 2005; Williams, Fuller, & Sridharan 2005), UC and compact 
H II regions (Chini, Kruegel, & Wargau 1987; Wood & Churchwell 1989a,b; Ball et al. 1996; 
Martin-Hernandez, van der Hulst, & Tielens 2003; Crowther & Conti 2003; Alvarez et al. 
2004), and masers (Testi et al. 1994; De Buizer, Pina, & Telesco 2000; De Buizer 2003a, 
2005a; Ellington 2006; Wu et al. 2006). Relatively few studies (e.g., Zoonematkermani et 
al. 1990; Becker et al. 1994; Giveon et al. 2005a; White, Becker, & Helfand 2005; Beltran 
et al. 2006) combined IR and radio observations to determine the characteristics of massive 
star- formation on a Galactic scale by analyzing large samples (> 100 sources). One of the 
major drawbacks of these larger samples was the poor angular resolution the IR surveys 
used - IRAS with ~ 1', and the Midcourse Space Experiment (MSX) with 18.3", compared 
to the resolution of the radio surveys (typically 1" to a few arcseconds). Not only does 
the resolution mismatch severely limit comparison between the IR and radio data, it also 
increases the risk of false positives in IR-radio samples. 

This work presents high-resolution diffraction-limited MIR observations of 18 UC H II 
region candidates. These objects are part of a sample of 687 candidates selected by matching 
IR and radio maps of the first Galactic quadrant (Giveon et al. 2005a). The 18 sources 
sample almost triples the number of long wavelength (> 18/im) MIR observations of UC H II 
regions observed at high-resolution (De Buizer et al. 2002a,b, 2003b, 2005b, and Kraemer 
et al. 2003). The IR maps for finding these 687 matches were taken from the MSX archive 
(Egan, Price, & Kraemer 2003), presenting a relative improvement in angular resolution 
compared to IRAS. The radio maps were generated from VLA observations (White, Becker, 
& Helfand 2005) at an unprecedented sensitivity (> 90% complete for sources with F 5G h z > 3 
mJy) and angular resolution (6"). This sample is the most complete of its kind so far, due 
to the quality and spatial extent of the radio maps. Giveon et al. (2005a) showed that the 
radio survey is > 90% complete in detecting all embedded O stars across the Galaxy in the 
survey's area, assuming ionization-bounded nebulae. The still relatively poor resolution of 
MSX can lead to confusion regarding IR-radio associations and prohibits the study of source 
morphologies. With higher IR resolution, it becomes feasible to study the relationships 
between the central stellar sources, the ionized gas, and the surrounding dust. The goal of 
the present observations is to verify our selection of UC H II regions by determining their 
sizes, luminosities, morphologies, and stellar types of the ionizing stars. 

The SPITZER program, Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE and GLIMPSE II), will provide maps of large parts of the Galactic plane < 
65°, |6| < 1°; Benjamin et al. 2003) at wavelengths < 8 fim. The major advantages of this 
survey compared to the MSX survey is its higher angular resolution (~ 2"), and its much 
higher sensitivity (0.4 mJy compared to 0.1 Jy at 8/zm). While GLIMPSE will be ideal in 
confirming the radio and IR coincidences, and will help improve our understanding of the 
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morphologies of our full sample of candidates, observations at longer wavelengths are needed 
for evaluating the dust properties. 

The new Keck observations present a factor of 50 improvement in angular resolution 
compared to MSX. While we have begun our study in a limited sub-sample, especially for 
our longest wavelength (24.5 /im), in the future, we hope to include more sources in the 
analysis presented in this paper. In §2 we describe the observations and present the sample. 
In §3 we present flux density maps, dust temperature maps, and dust optical depth maps of 
our sample, and we derive MIR luminosities and spectral types. In §4 we discuss correlations 
between the measured and calculated properties. In §5 we present our conclusions. 

2. Observations 

Observations were performed in two half-nights in August 2004 on the Keck I 10 meter 
telescope on Mauna Kea, using the Long- Wavelength Spectrometer (LWS). We obtained 
images at 11.7 /im (AA = 1.1 /jm), 17.65 /im (AA = 0.85 /mi), and 24.5 /im (AA = 0.74 
/im). The LWS detector is a 128 x 128 pixel Si: As array manufactured by Boeing. LWS has 
a field of view of 10.84" x 10.84", with a scale of 0.0847" pixel" 1 . Sky and optics background 
offsets were subtracted using secondary chopping throws of 10" 15" at 2.5 Hz and by nodding 
the telescope every 10 seconds. Frame times of 10-30 ms were used for all observations. The 
standard stars (3 Oph, a Aql, HD 2486, and HD 25477 were observed in the 11.7 /im and the 
17.65 /im bands, and in a range of airmasses, during the first night. The standard stars (5 
Oph and 7 Aql were observed in all 3 bands, and in a range of airmasses, during the second 
night. The FWHM of the standard stars was 0.3" and 0.4" at 11.7 and 17.65 /im, respectively, 
in the first night, and 0.4", 0.5", and 0.6" at 11.7, 17.65, and 24.5 /im, respectively, in the 
second night. These values are very close to the diffraction limits - 0.3", 0.4", and 0.6" at 
11.7, 17.65, and 24.5 /im, respectively. 

We designed our list of objects in a manner that allowed an efficient source acquisition, 
combining the sample properties and the observing methodology with the LWS instrument. 
The original sample from Giveon et al. (2005a) spans longitudes 340°-42° and latitudes 
|6| < 0.4°. We chose from the original sample all point sources (< 5" at 5 GHz) with IR- 
radio matching reliability > 95% (see definition in Giveon et al. 2005a) to minimize the 
probability for chance coincidences. We chose sources with rising IR spectrum (8-21 /im), to 
include only sources with dominant IR continuum associated with blackbody temperatures 
< 200 K. We limited ourselves to longitudes > 10° to overlap with the GLIMPSE I survey. 
We selected only sources that were isolated (within 30") in the MSX maps to avoid confusion 
of the chop-nod procedure, used to subtract the background from the images. Since time 
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allowed, two additional sources were added to the observations from a matching of the 
MSX6C catalog and a 1.4 GHz radio survey (Giveon et al. 2005b), which included later 
targets. 

These criteria resulted in a list of 54 targets at 5 GHz. From these, 22 fields were 
observed, out of which only 5 yielded no MIR source within ~ 5" of the radio positions. 
The two 1.4 GHz sources yielded one LWS detection, summing up to 24 observed fields with 
18 detections. All of these 24 sources have MSX-radio counterparts within < 12" and 16 
of them have their counterparts within < 4". We suspect that the non-detections result 
from the large difference in angular resolution between LWS and MSX (x50), combined 
with the small (10") field of view of LWS. This either caused missing the source, or a false 
match between the IR source and a radio source, or if the match is real, the IR source may 
be a low surface brightness source, below the detection limit of LWS. Because the targets 
were invisible in the optical, we offset the telescope to the target fields from visible stars 
with known coordinates, resulting in pointing accuracy typically better than 1". This makes 
the non-detections less likely to be a result of the telescope. The GLIMPSE survey will 
provide higher-resolution images of our entire survey of 687 sources, and will allow a better 
understanding of the relation between the IR and the radio emission from these regions. 
Table 1 lists the 24 RA-sorted radio positions (the centers of the LWS fields) in columns 
(2)-(3); the MIR total flux densities at the corresponding MSX bands (8.3, 12.1, 14.7, and 
21.3 /im) in columns (4)-(7); the integrated radio flux densities at 1.4 and 5 GHz (White, 
Becker, k Helfand 2005; Helfand, Becker, & White 2006) in columns (8)-(9); and an LWS 
detection flag in column (10). The sources from the 1.4 GHz survey have only 1.4 GHz data 
listed as their radio flux densities. 



3. Results and Data Reduction 

3.1. Flux Densities 

Standard stars were observed throughout the nights in all three bands and at varying 
airmass values. The best-fit photometric solution was obtained for all stars of a given night 
simultaneously. The color terms are altogether insignificant, and only small airmass cor- 
rections, which increase with wavelength, are required. The absolute photometric accuracy 
was derived by summing in quadrature the standard deviation of the standard stars flux 
densities from the least-squares airmass fit and the Poisson error of the measured counts. 
The 1(7 accuracies are typically 4% at 11.7 /im, 5% at 17.65 /im, and 20% at 24.5 pm. In 
some cases, the subtraction of the chop-nod images resulted in artificial offsets, making parts 
of the image have negative counts. These residual offsets were removed manually by setting 



Table I. Observed sources and their total MSX and radio flux densities 
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Note. — UC HII regions candidates observed, sorted by RA. Listed arc the radio positions, total MIR. flux densities of the IR. counterparts measured by 
the MSX point-source catalog (8.3, 12.1, 14.7, and 21.3 M m )> integrated radio flux densities (1.4 and 5 GHz), and LWS detection flag. 

a Morc reliable flux densities from a new multi-configuration VLA observations (Hclfand, Becker, &c White 2006J. Sec http://third.ucllnl.org/gps/indcx.html 
for details. 
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regions with no apparent emission to zero in each image. A cut-off of 3a above the sky level 
was applied to all images, and only local peaks above this cut-off are considered. Total flux 
densities were calculated by summing-up all pixels above the 3a cut-off and applying the 
absolute calibration. 

Figures 1, 2, and 3 show the LWS spectral energy distributions (SEDs) of our sources 
combined with GLIMPSE (flux densities at 3.6, 4.5, 5.8, and 8.0 /zm; see Table 2) and MSX 
data (Table 1): empty circles indicate GLIMPSE data points, filled circles indicate MSX data 
points, and triangles indicate LWS data points. These SEDs are not extinction corrected. 

Most SEDs demonstrate the consistency of our calibrated flux densities and those from 
GLIMPSE and MSX, mainly for sources that are MSX point sources and are confined to 
the field of view in the LWS images (see LWS images, Figures 4 - 31). MSX had a beam 
of 18.3", while LWS has a field of view of only 10" across, which leads to some inconsis- 
tencies for sources with emission extending the LWS field of view (e.g., 11.94545—0.03634 
in Figure 1). These plots show that our sources have rising MIR spectra typical of warm 
dust peaking at longer wavelengths (60 — 100 /im), but we note the irregularly flat SED of 
30.04343-0.14200 (Figure 2) between 10 and 21 /im, consistently seen in both the LWS and 
the MSX data. An extinction correction will raise the flux density at the shorter wavelengths, 
making 30.04343—0.14200 an inverted MIR spectrum source, indicative of a hotter source 
(see Table 3 in §3.2 for temperature information). This source is the weakest radio emitter 
in our sample, and it has a non- inverted radio spectrum (i.e., its 1.4 GHz flux density is 
greater than its 5 GHz flux density; see Table 1). This is possibly not an UC H II region but 
another type of source, such as an AGB star, a Wolf-Rayet star, or a supernova remnant. 
Additional observations are required to determine the nature of this source. 

Table 2 lists GLIMPSE counterparts of our 18 LWS sources. The GLIMPSE name 
designations are given in column (2); the flux densities (not corrected for extinction) at 3.6, 
4.5, 5.8, and 8.0 /iin in Janskys are listed in columns (3)-(6); and the distances in arcseconds 
from the LWS source are given in column (7). Table 3 lists the detected sources and their 
LWS flux densities: multiple peaks are numbered in column (2); peak flux densities in the 
three bands observed are given in columns (3), (6), and (9); peak locations (RA and Dec 
offsets from the coordinates of the field centers given in Table 1) are listed in columns (4), 

(7) , and (10); and their total flux densities in the bands observed are listed in columns (5), 

(8) , and (11). For fields with multiple sources, we label them '1', '2', '3', etc, in decreasing 
order of peak flux density at 11.7 fim. The same numerical labels at different frequencies 
correspond to the same peaks. 
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Fig. 1. — Combined GLIMPSE-MSX-LWS SEDs of our sources with their la uncertainties: 
empty circles indicate GLIMPSE data points, filled circles indicate MSX data points, and 
triangles indicate LWS data points. Source identification is given in each panel. 



- 9- 




10 15 

X [|im] 



10 15 

X [|im] 



Fig. 2. — Combined GLIMPSE-MSX-LWS SEDs of our sources with their la uncertainties: 
empty circles indicate GLIMPSE data points, filled circles indicate MSX data points, and 
triangles indicate LWS data points. Source identification is given in each panel. 
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Fig. 3. — Combined GLIMPSE-MSX-LWS SEDs of our sources with their la uncertainties: 
empty circles indicate GLIMPSE data points, filled circles indicate MSX data points, and 
triangles indicate LWS data points. Source identification is given in each panel. 
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Table 2. GLIMPSE flux densities of our detected sources 
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Note. — GLIMPSE counterparts of our observed sources. We list the GLIMPSE name designation, the 3.6, 4.5, 
5.8, and 8.0 fim flux densities in Janskys and the distance D in arcseconds from the LWS source. 
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Flux density maps in the observed bands are shown in Figures 4-31 with sub-structures 
indicated by a cross as listed in Table 3. The images shown are smoothed to make their 
resolution compatible with each other for the subsequent analysis (§3.2). The smoothed 
resolution is given in the individual images. All images were rotated, so that north is up 
and east is to the left. The directions of increasing Galactic longitude I and latitude b are 
indicated in each map. Eight of our sources have radio maps with similar resolution in the 
literature: The UC H II regions 11.11198-0.39795, 28.28875-0.36359, 37.87411-0.39866, 
and 111.28293-0.66355 (Kurtz, Churchwell, & Wood 1994) selected using IRAS colors; the 
radio-selected sources 25.39918-0.14081, 33.91585+0.11111, and 35.46832+0.13984 (Fey, 
Spangler, & Cordes 1991); and 27.18725-0.08095 (Sridharan et al. 2002) selected as a 
high-mass protostellar object using FIR, radio continuum, and molecular line data. The 
radio maps show a good overall match in shape when compared to our LWS images. 



3.2. Dust Temperatures 

Dust color temperatures maps for each source were generated from the flux density 
maps. When a source is observed in three bands, both Ti 7.65/11. 7 and T24.5/17.Q5 are calcu- 
lated. These temperatures are usually not the same, because the sources are not necessarily 
perfect blackbodies. Detailed modeling of the dust properties, the geometry, and the pos- 
sible spectral features in the bands observed for each source are beyond the scope of this 
paper. Even though our calculation is less accurate, the color temperatures and the optical 
depths (§3.3) are used to understand the spatial trends in our sources, and their face values 
are secondary. See Li & Chen (1996), De Buizer et al. (2002a), and De Buizer (2005a) for 
similar analysis. 

A few preparatory steps are required in order to construct the temperature maps prop- 
erly. First, all images were aligned relative to the 11.7 jum images by employing a centroid- 
fitting algorithm to the peak or peaks. The images were then smoothed to the resolution of 
the image with the poorest resolution in each pair of flux density maps (17.65 fim or 24.5 
/im). The smoothing was done by convolving the images with the point-spread function 
(PSF) of one of the reference stars in the corresponding band. For each of the two tem- 
peratures, we have also convolved the poorer resolution image with the PSF of the better 
resolution band. For example, the 11.7 /zm image was convolved with the 17.65 /zm PSF, and 
the 17.65 /im image was convolved with the 11.7 /xm PSF. This ensured that both images 
had the same effective resolution (De Buizer et al. 2002a). A 3c cut-off above the noise level 
was applied to all images, to ensure using only pixels with minimally good signal-to-noise. 
Temperatures were then calculated only for diffraction limited spatial resolution elements 
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which were above this cut-off at both bands. 

We have tested the robustness of our temperature maps by shifting one of the flux density 
maps of each temperature by ±3 pixels (0.25") in RA and Dec - an unlikely large uncertainty 
in the image alignment. The absolute values of the temperature maxima changed by as much 
as 30%, but the main components remained unchanged in their shapes and relative positions. 
The average temperatures are much less sensitive to misalignments, and for the above shifts 
in the flux density maps, they typically vary by < 1%. This test also shows that aligning 
the images using the peak intensity does not bias the temperature and optical depth peaks. 

Temperatures were computed pixel by pixel using a look-up table containing the calcu- 
lated flux density ratios as a function of temperatures, assuming a modified blackbody with 
the dust emissivity law of Draine (2003) for Ry = 3.1 2 , not including foreground extinction, 
and taking into account the filter transmission and the atmospheric transmission through the 
bandpass. According to the Draine (2003) emissivity law, the corresponding optical depths 
at the central wavelengths of our bands (incorporating the filter bandpasses) have the rela- 
tions Ti7. 65 /tii.7 = 0.71 and T 2 4.5/r 1765 = 0.64. The typical temperature uncertainty from 
the look-up table is 1-2 K for T 17 . 65 / 1L7 , and 1-3 K for T 24 . 5 /i 7 . 6 5. Table 3 lists the labels of 
multiple temperature peaks in column (2); the values of the temperature peaks in columns 
(3) and (6); the locations of the peaks (RA and Dec offsets from the coordinates of the field 
centers given in Table 1) in columns (4) and (7); and the mean temperatures in columns 
(5) and (8) for each source. Multiple temperature peaks are labeled by the upper-case let- 
ters 'A', 'B', £ C, etc (column 2), and are sorted by decreasing T 17 Q 5 / n 7 . The upper-case 
alphabetical labels at different frequencies correspond to the same peaks. We note again the 
higher temperature of the exceptional source 30.04343—0.14200 reflecting the unusual flat 
SED of this source. We suspect that this source is not an UC HII region (see MIR SED in 
Figure 2, §3.1). Figures 4-31 show the temperature maps as contour plots overlaid on the 
flux density maps for all of our sources. Peak temperatures are indicated by a cross, with 
sub-peaks indicated as in Table 3, where 'A' is the highest peak. 



2 The extinction curve is available at http://www.astro.princeton.edu/~draine/dust/dustmix.html 
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Table 3. Temperature Maps Data 





Peak 


17 fi^ / 1 1 1 
1 1 .\>D / 11. 1 


Offset (Aa,A<5) 


< ^17.65/11.7 > 


T-imax 

OA t;/l7 fix. 

/LI. Oil 


Offset (Aq,A5) 


< ^24.5/17.65 > 


Source 


Label 


(K) 


(arcsec) 


(K) 


(K) 


(arcsec) 


(K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


11.11198-0.39795 




184 


(+3.7,-3.4) 


143 


81 


(+0.6,-0.7) 


72 


11.94545-0.03634 


A 








81 


(+0.7,+1.9) 


71 




B 


159 


(-1.3.+1.7) 


160 


78 


(-0.9.+2.5) 






C 








75 


(-1.4,-0.2) 




18.71179+0.00085 




154 


(+0.4.+0.6) 


149 








19.75611-0.12775 




159 


(+0.5,-0.1) 


142 


91 


(+0.6,-0.2) 


77 


21.38654-0.25346 




159 


(-0.4,+0.5) 


152 








25.39918-0.14081 


A 


158 


(-2.1,-0.2) 


128 










B 








101 


(+0.7,0.0) 


76 




C 








98 


(+0.7,-1.2) 






D 








96 


(+0.9,-3.9) 




25.80211-0.15640 




143 


(-0.8.+0.6) 


128 


101 


(-1.4.+1.2) 


85 


27.18725-0.08095 


A 


140 


(+0.5,+0.2) 


135 


62 


(0.0, +0.5) 


64 




B 








84 


(+1.3,+3.2) 






C 








63 


(+0.4,-0.9) 




28.28875-0.36359 


A 


174 


(+1.0,+1.3) 


119 


104 


(+1.1, +0.8) 


87 




B 








95 


(-2.5,-3.3) 




30.04343-0.14200 


A 


327 


(-0.3.+0.4) 


251 










B 


243 


(-0.2,-1.8) 










30.86744+0.11493 


A 


215 


(-0.3.+0.3) 


183 










B 


212 


(-0.6,-1.4) 










30.58991-0.04231 




229 


(+0.3,-0.1) 


213 








30.66808-0.33134 




160 


(+0.8,0.0) 


147 








33.91585+0.11111 


A 


160 


(+0.6,+0.6) 


138 


109 


(+0.3,+0.3) 


87 




B 








100 


(+1.6,-1.9) 






C 








99 


(+0.6,-1.4) 




33.81104-0.18582 


A 


231 


(0.0,-0.5) 


150 


89 


(0.0,-0.8) 


70 




B 


138 


(+0.3,-2.1) 




80 


(+0.4,-1.9) 




35.46832+0.13984 




143 


(-1.1, +0.7) 


131 








37.87411-0.39866 


A 


195 


(0.0,-0.2) 


149 


85 


(+0.5,+0.7) 


78 




B 


164 


(+0.7,-2.8) 












C 


162 


(-3.3,-1.6) 












D 








88 


(-1.9.+1.5) 




111.28293-0.66355 


A 


133 


(+2.2,+2.6) 


130 










B 


117 


(-0.2.+3.0) 











Note. — Temperature peaks, positions of the peaks, and temperature means for our sources. The listed values are from our best alignment 
of the flux density maps. The typical lcr uncertainties due to a possible misalignment arc 30% for the peak temperatures and 1% for the 
mean temperatures. Some of the fields contain multiple temperature peaks, which are labeled 'A', 'B', 'C, etc., and are indicated in the 
maps in Figures 4-31. 
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6 4 2 -2 -4 4 2 2 -4 -6 
RA offset (a rosec) R A offset (arcsec) 

Fig. 4. — Maps of the source 11.11198—0.39795: (a) The 11.7 fim flux density image 
smoothed to a resolution of 0.4" with overplotted contours. The directions of increasing 
Galactic longitude I and latitude b are indicated, (b) The 17.65 /im flux density image with 
overplotted contours, (c) Tns5/n.7 plotted as contours overlaid on the 17.65 /xm flux density 
map. The region near the flux density peak is the temperature minimum, (d) rn. 7 plotted as 
contours overlaid on the 17.65 /im flux density map. Sharp edges result from array rotation 
to orient north up and east to the left. Flux density contour levels are 90, 75, 50, 25, 10, 5, 
and 1 percent of the map peak - 1.12, 0.93, 0.62, 0.31, 0.12, 0.06, 0.01 Jansky/arcsec 2 at 11.7 
/im with a peak of 1.24 Jansky/arcsec 2 and 5.1, 4.3, 2.8, 1.4, 0.6, 0.3, 0.06 Jansky/arcsec 2 
at 17.65 /im with a peak of 5.7 Jansky/arcsec 2 . Temperature contour levels are 95, 90, 80, 
70, 60, and 50 percent of the map peak (184 K) - 175, 166, 147, 129, 110, 92 K. Optical 
depth contour levels are 75, 50, 25, 10, 5, and 1 percent of the map peak (0.07) - 0.053, 
0.035, 0.018, 0.007, 0.004, 0.0007. Peaks are indicated by a cross. Multiple peaks are sorted 
numerically or alphabetically (see Tables 3 and 4). 
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6 4 2 -2 -4 4 2 -2 -4 -6 
RA offset (arcsec) RA offset (arosec) 

Fig. 5.— Maps of the source 11.11198-0.39795: (a) The 17.65 flux density image 
smoothed to a resolution of 0.6" with overplotted contours. The directions of increasing 
Galactic longitude I and latitude b are indicated, (b) The 24.5 fim flux density image 
with overplotted contours, (c) T24.5/17.Q5 plotted as contours overlaid on the 24.5 /im flux 
density map. (d) r 2 4.5 plotted as contours overlaid on the 24.5 fim flux density map. Sharp 
edges result from array rotation to orient north up and east to the left. The 17.65/zm flux 
density contour levels are 5.1, 4.3, 2.8, 1.4, 0.6, 0.3, 0.06 Jansky/arcsec 2 with a peak of 5.7 
Jansky/arcsec 2 . The 24.5/zm flux density contour levels are 19.8, 16.5, 11.0, 5.5, 2.2, 1.1, 0.2 
Jansky/arcsec 2 with a peak of 22 Jansky/arcsec 2 . Temperature contour levels are 77, 73, 65, 
57, 49, 41 K, with a temperature peak of 81 K. Optical depth contour levels are 0.53, 0.35, 
0.18, 0.07, 0.04, 0.007 with a peak of 0.7. Peaks are indicated by a cross. Multiple peaks are 
sorted numerically or alphabetically (see Tables 3 and 4). 
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Fig. 6. — Maps of the source 11.94545—0.03634: (a), (b), (c), and (d) are the same as in 
Figure 4. The 11.7/xm flux density contour levels are 1.11, 0.92, 0.62, 0.31, 0.12, 0.06, 0.01 
Jansky/arcsec 2 with a peak of 1.23 Jansky /arcsec 2 . The 17.65/im flux density contour levels 
are 2.36, 1.96, 1.31, 0.65, 0.26, 0.13, 0.03 Jansky/arcsec 2 with a peak of 2.62 Jansky/arcsec 2 . 
Temperature contour levels are 151, 143, 127, 111, 95, 80 K, with a temperature peak of 
159 K. Optical depth contour levels are 0.0049, 0.0033, 0.0016, 0.0007, 0.00033, 6.5 • 10~ 5 
with a peak of 0.0065. Sharp edges result from array rotation to orient north up and east to 
the left. In the temperature map, the region near the flux density peak is the temperature 
minimum. 
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Fig. 7. — Maps of the source 11.94545—0.03634: (a), (b), (c), and (d) are the same as in 
Figure 5. The 17.65/zm flux density contour levels are 2.36, 1.96, 1.31, 0.65, 0.26, 0.13, 0.03 
Jansky/arcsec 2 with a peak of 2.62 Jansky/arcsec 2 . The 24.5/zm flux density contour levels 
are 8.93, 7.44, 4.96, 2.48, 0.99, 0.50, 0.10 Jansky/arcsec 2 with a peak of 9.92 Jansky/arcsec 2 . 
Temperature contour levels are 77, 73, 65, 57, 49, 41 K, with a temperature peak of 81 K. 
Optical depth contour levels are 0.25, 0.17, 0.08, 0.03, 0.02, 0.003 with a peak of 0.33. Sharp 
edges result from array rotation to orient north up and east to the left. 
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Fig. 8. — Maps of the source 18.71179+0.00085: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/im flux density contour levels are 0.53, 0.44, 0.29, 0.15, 0.06, 0.03, 
0.006 Jansky/arcsec 2 with a peak of 0.59 Jansky/arcsec 2 . The 17.65//m flux density contour 
levels are 1.7, 1.4, 0.9, 0.5, 0.2, 0.1, 0.02 Jansky/arcsec 2 with a peak of 1.89 Jansky/arcsec 2 . 
Temperature contour levels are 146, 137, 123, 108, 92, 77 K, with a temperature peak of 154 
K. Optical depth contour levels are 0.0023, 0.0015, 0.0008, 0.0003, 0.0002, 3.0 • 10~ 5 with a 
peak of 0.003. 
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Fig. 9. — Maps of the source 19.75611—0.12775: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/zm flux density contour levels are 5.0, 4.2, 2.8, 1.4, 0.6, 0.3, 0.06 
Jansky/arcsec 2 with a peak of 5.6 Jansky /arcsec 2 . The 17.65/iin flux density contour levels 
are 16.6, 13.8, 9.2, 4.6, 1.8, 0.9, 0.2 Jansky/arcsec 2 with a peak of 18.4 Jansky/arcsec 2 . 
Temperature contour levels are 151, 143, 127, 111, 95, 80 K, with a temperature peak of 159 
K. Optical depth contour levels are 0.022, 0.015, 0.007, 0.003, 0.001, 0.0003 with a peak of 
0.029. 
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Fig. 10.— Maps of the source 19.75611-0.12775: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65//m flux density contour levels are 16.6, 13.8, 9.2, 4.6, 1.8, 0.9, 0.2 
Jansky/arcsec 2 with a peak of 18.4 Jansky/arcsec 2 . The 24.5/zm flux density contour levels 
are 42.3, 35.3, 23.5, 11.8, 4.7, 2.4, 0.5 Jansky/arcsec 2 with a peak of 47 Jansky/arcsec 2 . 
Temperature contour levels are 86, 82, 73, 64, 55, 46 K, with a temperature peak of 91 K. 
Optical depth contour levels are 0.45, 0.30, 0.15, 0.06, 0.03, 0.006 with a peak of 0.60. 
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Fig. 11. — Maps of the source 21.38654—0.25346: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/zm flux density contour levels are 1.22, 1.01, 0.68, 0.34, 0.14, 0.07 
Jansky/arcsec 2 with a peak of 1.35 Jansky /arcsec 2 . The 17.65/im flux density contour levels 
are 3.6, 3.0, 2.0, 1.0, 0.4, 0.2 Jansky/arcsec 2 with a peak of 4 Jansky/arcsec 2 . Temperature 
contour levels are 151, 143, 127 K, with a temperature peak of 159 K. Optical depth contour 
levels are 0.0082, 0.0055, 0.0027, 0.0011, 0.0005 with a peak of 0.011. 
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Fig. 12.— Maps of the source 25.39918-0.14081: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/xm flux density contour levels are 0.9, 0.75, 0.5, 0.25, 0.1, 0.05, 
0.01 Jansky/arcsec 2 with a peak of 1.0 Jansky/arcsec 2 . The 17.65/im flux density contour 
levels are 3.8, 3.1, 2.1, 1.1, 0.4, 0.2, 0.04 Jansky/arcsec 2 with a peak of 4.2 Jansky/arcsec 2 . 
Temperature contour levels are 150, 142, 126, 111, 95, 79 K, with a temperature peak of 158 
K. Optical depth contour levels are 0.014, 0.009, 0.005, 0.002, 0.001, 0.0002 with a peak of 
0.019. Sharp edges result from array rotation to orient north up and east to the left. 
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Fig. 13. — Maps of the source 25.39918—0.14081: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65/xm flux density contour levels are 3.8, 3.1, 2.1, 1.1, 0.4, 0.2, 0.04 
Jansky/arcsec 2 with a peak of 4.2 Jansky/arcsec 2 . The 24.5/zm flux density contour lev- 
els are 12.6, 10.5, 7.0, 3.5, 1.4, 0.7, 0.1 Jansky/arcsec 2 with a peak of 14 Jansky/arcsec 2 . 
Temperature contour levels are 96, 91, 81, 71, 61, 51 K, with a temperature peak of 101 K. 
Optical depth contour levels are 0.60, 0.40, 0.20, 0.08, 0.04, 0.008 with a peak of 0.80. Sharp 
edges result from array rotation to orient north up and east to the left. 
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Fig. 14. — Maps of the source 25.80211—0.15640: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/zm flux density contour levels are 3.7, 3.1, 2.1, 1.0, 0.4, 0.2, 0.04 
Jansky/arcsec 2 with a peak of 4.1 Jansky /arcsec 2 . The 17.65//m flux density contour levels 
are 14.3, 11.9, 8.0, 4.0, 1.6, 0.8, 0.2 Jansky/arcsec 2 with a peak of 15.9 Jansky/arcsec 2 . 
Temperature contour levels are 136, 129, 114, 100, 86, 72 K, with a temperature peak of 143 
K. Optical depth contour levels are 0.029, 0.020, 0.010, 0.004, 0.002, 0.0004 with a peak of 
0.039. 
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Fig. 15. — Maps of the source 25.80211—0.15640: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65/iin flux density contour levels are 14.3, 11.9, 8.0, 4.0, 1.6, 0.8, 0.2 
Jansky/arcsec 2 with a peak of 15.9 Jansky/arcsec 2 . The 24.5/zm flux density contour levels 
are 36.0, 30.0, 20.0, 10.0, 4.0, 2.0, 0.4 Jansky/arcsec 2 with a peak of 40 Jansky/arcsec 2 . 
Temperature contour levels are 96, 91, 81, 71, 61, 51 K, with a temperature peak of 101 K. 
Optical depth contour levels are 0.26, 0.18, 0.09, 0.04, 0.02, 0.003 with a peak of 0.35. 
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Fig. 16. — Maps of the source 27.18725—0.08095: (a), (b), (c), and (d) are the same as in 
Figure 4. The 11.7/iin flux density contour levels are 0.39, 0.32, 0.22, 0.11, 0.04, 0.02, 0.004 
Jansky/arcsec 2 with a peak of 0.43 Jansky /arcsec 2 . The 17.65/zm flux density contour levels 
are 1.75, 1.46, 0.97, 0.49, 0.19, 0.10, 0.02 Jansky/arcsec 2 with a peak of 1.94 Jansky/arcsec 2 . 
Temperature contour levels are 133, 126, 112, 98, 84, 70 K, with a temperature peak of 140 
K. Optical depth contour levels are 0.0038, 0.0025, 0.0013, 0.0005, 0.0003, 5 • 10~ 5 with a 
peak of 0.005. 
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Fig. 17.— Maps of the source 27.18725-0.08095: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65//m flux density contour levels are 1.75, 1.46, 0.97, 0.49, 0.19, 0.10, 
0.02 Jansky/arcsec 2 with a peak of 1.94 Jansky/arcsec 2 . The 24.5/zm flux density contour 
levels are 16.2, 13.5, 9.0, 4.5, 1.8, 0.9, 0.2 Jansky/arcsec 2 with a peak of 18 Jansky/arcsec 2 . 
Temperature contour levels are 80, 76, 67, 59, 50, 42 K, with a temperature peak of 84 K. 
Optical depth contour levels are 5.3, 3.5, 1.8, 0.7, 0.35, 0.01 with a peak of 7. 
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Fig. 18.— Maps of the source 28.28875-0.36359: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/zm flux density contour levels are 9.2, 7.6, 5.1, 2.6, 1.0, 0.5, 0.1 
Jansky/arcsec 2 with a peak of 10.2 Jansky /arcsec 2 . The 17.65/im flux density contour levels 
are 23.4, 19.5, 13.0, 6.5, 2.6, 1.3, 0.3 Jansky/arcsec 2 with a peak of 26 Jansky/arcsec 2 . 
Temperature contour levels are 133, 126, 112, 98, 84, 70 K, with a temperature peak of 174 
K. Optical depth contour levels are 0.09, 0.06, 0.03, 0.012, 0.006, 0.001 with a peak of 0.12. 
Sharp edges result from array rotation to orient north up and east to the left. The central 
bulb between the two peaks in the r n 7 map is a local minimum. 
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Fig. 19. — Maps of the source 28.28875—0.36359: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65/iin flux density contour levels are 23.4, 19.5, 13.0, 6.5, 2.6, 1.3, 0.3 
Jansky/arcsec 2 with a peak of 26 Jansky /arcsec 2 . The 24.5/zm flux density contour levels 
are 54.0, 45.0, 30.0, 15.0, 6.0, 3.0, 0.6 Jansky/arcsec 2 with a peak of 60 Jansky/arcsec 2 . 
Temperature contour levels are 99, 94, 83, 73, 62, 52 K, with a temperature peak of 104 K. 
Optical depth contour levels are 0.43, 0.29, 0.14, 0.06, 0.03, 0.006 with a peak of 0.57. Sharp 
edges result from array rotation to orient north up and east to the left. 
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3.3. Emission Optical Depths 



The measured MIR intensity emitted from dust grains is related to their temperature 



where B u (Td) is the Planck function at dust temperature To, t v is the emission optical 
depth of the dust, and (1 — e~ Tv ) is called the emissivity function. This emission is typically 
optically thin in UC H II regions (see e.g., De Buizer et al. 2002a; and De Buizer 2005a), 
which are small, hot sources, so the emissivity function may be reduced to simply r v . The 
optical depth may be offset by a proportion constant due to absorption from a much larger 
and cooler cloud. We assume that the absorption is uniform over our fields of view, since it 
arises on scales much larger than our sources. The optically thinness assumption does not 
hold in the case of 27.18725—0.08095, with mean T24.5 > 1. The peak T24.5 of the source 
37.87411—0.39866 is also larger than 1, but the average T24.5 is still < I. 

We have generated emission optical depth maps at 11.7 /im and at 24.5 /im, pixel 
by pixel, using their corresponding flux density and dust temperature maps. For sources 
observed only in two bands, only 711.7 wa s calculated. Flux densities were translated into 
intensities by normalizing them by their angular pixel size. We have tested the robustness of 
our optical depth maps by shifting one of the flux density maps by ±3 pixels in RA and Dec, 
as was done when creating the temperature maps. The absolute values of the optical depth 
maxima changed by as much as x2, but the main components remained unchanged in their 
shapes and relative positions. The average optical depths are less sensitive to misalignments, 
and for the above shifts in the flux density maps, they typically vary by 50%. 

Table 4 lists the calculated 11.7 /im and 24.5 fim optical depth data: labels of multiple 
peaks in column (2); the values of the peaks in columns (3) and (6); the locations of the peaks 
(RA and Dec offsets from the coordinates of the field centers given in Table 1) in columns (4) 
and (7); and mean optical depth in columns (5) and (8). Some of the fields contain multiple 
optical depth peaks, which are labeled by the lower-case letters 'a', 'b', 'c', etc, sorted by 
decreasing r n 7. The same lower-case alphabetical labels at different frequencies correspond 
to the same peaks. Figures 4-31 show optical depth maps as contour plots overlaid on 
the 17.65 /xm or 24.5 fim flux density maps for all of our sources. Peak optical depths are 
indicated by a cross, with sub-peaks indicated as in Table 4, where 'a' is the highest peak. 

In general, T24.5 > Tn.7 for all of our sources. This is because first, at 24.5 fim we can 
see deeper into the region compared to 11.7 fim. Second, 24.5 fim probes cooler dust, which 
is usually distributed further out from the central source, whereas 11.7 fim probes hotter 
dust further in (T17.65/11.7 > T24.5717.65 as m Table 3). The overall effect is that the 24.5 fim 
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Fig. 20. — Maps of the source 30.04343—0.14200: (a), (b), (c), and (d) are the same as in 
Figure 4. The 11.7 /im flux density image is smoothed to a resolution of 0.5". The 11.7/im 
flux density contour levels are 14.8, 12.3, 8.2, 4.1, 1.6, 0.8, 0.2 Jansky /arcsec 2 with a peak 
of 16.4 Jansky /arcsec 2 . The 17.65/xm flux density contour levels are 13.1, 11.0, 7.3, 3.7, 1.5, 
0.7, 0.1 Jansky /arcsec 2 with a peak of 14.6 Jansky/arcsec 2 . Temperature contour levels are 
311, 294, 262, 229, 196, 164 K, with a temperature peak of 372 K. Optical depth contour 
levels are 0.00113, 0.00075, 0.00038, 0.00015, 7.5 • 10~ 5 , 1.5 • 10~ 5 with a peak of 0.0015. 
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Fig. 21. — Maps of the source 30.86744+0.11493: (a), (b), (c), and (d) are the same as in 
Figure 4. The 11.7 /im flux density image is smoothed to a resolution of 0.5". The 11.7/im 
flux density contour levels are 6.3, 5.3, 3.5, 1.8, 0.7, 0.4, 0.07 Jansky/arcsec 2 with a peak of 
7.0 Jansky/arcsec 2 . The 17.65/zm flux density contour levels are 9.5, 7.9, 5.3, 2.6, 1.1, 0.5, 
0.1 Jansky/arcsec 2 with a peak of 10.5 Jansky/arcsec 2 . Temperature contour levels are 204, 
193, 172, 150, 129, 108 K, with a temperature peak of 215 K. Optical depth contour levels 
are 0.0026, 0.0017, 0.0009, 0.0003, 0.0002, 3.4 • 10~ 5 with a peak of 0.0034. 
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Fig. 22. — Maps of the source 30.58991—0.04231: (a), (b), (c), and (d) are the same as in 
Figure 4. The 11.7/iin flux density contour levels are 0.52, 0.43, 0.29, 0.15, 0.06, 0.03, 0.006 
Jansky/arcsec 2 with a peak of 0.58 Jansky /arcsec 2 . The 17.65/xm flux density contour levels 
are 0.78, 0.65, 0.44, 0.22, 0.09, 0.04, 0.009 Jansky/arcsec 2 with a peak of 0.87 Jansky/arcsec 2 . 
Temperature contour levels are 218, 206, 183, 160, 137, 114 K, with a temperature peak of 
229 K. Optical depth contour levels are 0.00021, 0.00014, 7 • 10" 5 , 2.8 • 10~ 5 , 1.4 • 10~ 5 , 
2.8 • 10~ 6 with a peak of 0.00028. 
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Fig. 23. — Maps of the source 30.66808—0.33134: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7//m flux density contour levels are 1.71, 1.43, 0.96, 0.48, 0.19, 0.10, 
0.02 Jansky /arcsec 2 with a peak of 1.91 Jansky/arcsec 2 . The 17.65/zm flux density contour 
levels are 5.8, 4.8, 3.2, 1.6, 0.6, 0.3, 0.06 Jansky/arcsec 2 with a peak of 6.4 Jansky/arcsec 2 . 
Temperature contour levels are 152, 144, 128, 112, 96, 80 K, with a temperature peak of 
160 K. Optical depth contour levels are 0.0073, 0.0049, 0.0024, 0.0010, 0.0005, 0.0001 with 
a peak of 0.0097. 
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Fig. 24. — Maps of the source 33.91585+0.11111: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/im flux density contour levels are 1.2, 1.0, 0.7, 0.3, 0.1, 0.07, 0.01 
J ansky/ arcsec 2 with a peak of 1.3 Jansky /arcsec 2 . The 17.65/zm flux density contour levels 
are 3.7, 3.1, 2.1, 1.0, 0.4, 0.2, 0.04 Jansky/arcsec 2 with a peak of 4.1 Jansky/arcsec 2 . Tem- 
perature contour levels are 152, 144, 128, 112, 96, 80 K, with a temperature peak of 160 K. 
Optical depth contour levels are 0.0089, 0.0060, 0.0030, 0.0012, 0.0006, 0.0001 with a peak 
of 0.012. Sharp edges result from array rotation to orient north up and east to the left. 
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Fig. 25. — Maps of the source 33.91585+0.11111: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65//m flux density contour levels are 3.7, 3.1, 2.1, 1.0, 0.4, 0.2, 0.04 
Jansky/arcsec 2 with a peak of 4.1 Jansky/arcsec 2 . The 24.5/xm flux density contour levels are 
8.1, 6.8, 4.5, 2.3, 0.9, 0.5, 0.09 Jansky/arcsec 2 with a peak of 9 Jansky/arcsec 2 . Temperature 
contour levels are 104, 98, 87, 76, 65, 54 K, with a temperature peak of 109 K. Optical depth 
contour levels are 0.060, 0.040, 0.020, 0.008, 0.004, 0.0008 with a peak of 0.08. Sharp edges 
result from array rotation to orient north up and east to the left. 
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Fig. 26. — Maps of the source 33.81104—0.18582: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/zm flux density contour levels are 7.6, 6.3, 4.2, 2.1, 0.8, 0.4, 0.08 
Jansky/arcsec 2 with a peak of 8.4 Jansky /arcsec 2 . The 17.65/zm flux density contour levels 
are 9.7, 8.1, 5.4, 2.7, 1.1, 0.5, 0.1 Jansky/arcsec 2 with a peak of 10.8 Jansky/arcsec 2 . Tem- 
perature contour levels are 219, 208, 185, 162, 139, 116 K, with a temperature peak of 231 
K. Optical depth contour levels are 0.01065, 0.0071, 0.0035, 0.0014, 0.0007, 0.0001 with a 
peak of 0.014. 
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Fig. 27.— Maps of the source 33.81104-0.18582: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65/zm flux density contour levels are 9.7, 8.1, 5.4, 2.7, 1.1, 0.5, 0.1 
Jansky/arcsec 2 with a peak of 10.8 Jansky/arcsec 2 . The 24.5/xm flux density contour levels 
are 24.3, 20.3, 13.5, 6.8, 2.7, 1.4, 0.3 Jansky/arcsec 2 with a peak of 27 Jansky/arcsec 2 . 
Temperature contour levels are 85, 80, 71, 62, 53, 44 K, with a temperature peak of 89 K. 
Optical depth contour levels are 0.40, 0.27, 0.13, 0.05, 0.03, 0.005 with a peak of 0.53. 
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Fig. 28.— Maps of the source 35.46832+0.13984: (a), (b), (c), and (d) are the same as in 
Figure 4. The 11.7 /xm flux density image is smoothed to a resolution of 0.5". The 11.7/xm 
flux density contour levels are 1.2, 1.0, 0.7, 0.3, 0.1, 0.07, 0.01 Jansky /arcsec 2 with a peak of 
1.38 Jansky /arcsec 2 . The 17.65/zm flux density contour levels are 5.6, 4.6, 3.1, 1.6, 0.6, 0.3, 
0.06 Jansky/arcsec 2 with a peak of 6.2 Jansky/arcsec 2 . Temperature contour levels are 136, 
129, 114, 100, 86, 72 K, with a temperature peak of 143 K. Optical depth contour levels are 
0.029, 0.019, 0.009, 0.004, 0.002, 0.0004 with a peak of 0.038. Sharp edges result from array 
rotation to orient north up and east to the left. 
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Fig. 29. — Maps of the source 37.87411—0.39866: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7/zm flux density contour levels are 7.4, 6.1, 4.1, 2.1, 0.8, 0.4, 0.08 
Jansky/arcsec 2 with a peak of 8.2 Jansky /arcsec 2 . The 17.65/zm flux density contour levels 
are 13.4, 11.2, 7.4, 3.7, 1.5, 0.7, 0.1 Jansky/arcsec 2 with a peak of 14.9 Jansky/arcsec 2 . 
Temperature contour levels are 185, 175, 156, 136 K, with a temperature peak of 195 K. 
Optical depth contour levels are 0.0056, 0.0038, 0.0019, 0.0008 with a peak of 0.0075. Sharp 
edges result from array rotation to orient north up and east to the left. 
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Fig. 30. — Maps of the source 37.87411—0.39866: (a), (b), (c), and (d) are the same as 
in Figure 5. The 17.65/im flux density contour levels are 13.4, 11.2, 7.4, 3.7, 1.5, 0.7, 0.1 
Jansky/arcsec 2 with a peak of 14.9 Jansky/arcsec 2 . The 24.5/zm flux density contour levels 
are 44.1, 36.8, 24.5, 12.3, 4.9, 2.5 0.5 Jansky/arcsec 2 with a peak of 49 Jansky/arcsec 2 . 
Temperature contour levels are 84, 79, 70, 62, 53, 44 K, with a temperature peak of 88 K. 
Optical depth contour levels are 0.90, 0.60, 0.30, 0.12, 0.06, 0.01 with a peak of 1.2. Sharp 
edges result from array rotation to orient north up and east to the left. 
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Fig. 31. — Maps of the source 111.28293—0.66355: (a), (b), (c), and (d) are the same as 
in Figure 4. The 11.7//m flux density contour levels are 0.74, 0.61, 0.41, 0.21, 0.08, 0.04, 
0.008 Jansky/arcsec 2 with a peak of 0.82 Jansky/arcsec 2 . The 17.65/xm flux density contour 
levels are 4.3, 3.6, 2.4, 1.2, 0.5, 0.2, 0.05 Jansky/arcsec 2 with a peak of 4.8 Jansky/arcsec 2 . 
Temperature contour levels are 126, 120, 106, 93, 80, 66 K, with a temperature peak of 133 
K. Optical depth contour levels are 0.017, 0.012, 0.006, 0.002, 0.001, 0.0002 with a peak of 
0.023. Sharp edges result from array rotation to orient north up and east to the left. 
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emission comes from a larger volume. We note that silicate absorption or emission from the 
strong features at 9.7 and 18 /im cannot be the reason for this effect in the optical depths. 
We have artificially changed the 11.7 and 17.65 /im flux densities to see the effect on the 
dust temperatures and optical depths. These tests show that the 18 /xm feature must be a 
stronger absorber or a weaker emitter than the part of the 9.7 /im feature affecting our data, 
in order not to reduce Tn.7 even more, which is unlikely. However, The 18 /im feature is 
always shallower in absorption than the 9.7 /im feature, and 11.7 /im is still well within the 
influence of the 9.7 /im feature. Also, UC H II regions generally have silicate absorption, in 
various depths from very shallow to very deep, but are hardly ever seen in emission. 

The optical depth peaks and the temperature peaks are often at different locations 
- the optical depth peaks are usually more closely-linked to the flux density peaks. The 
temperature peak is usually offset from the flux density peak. This means that having a 
large MIR flux density at a given location is usually due to a larger column of matter and/or 
higher density along that line-of-sight, and not necessarily due to the central ionizing source 
of the UC H II region. Such offsets have been predicted by blister-type models of H II regions, 
where the exciting star lies near the edge/outside the molecular cloud (Icke, Gatley, & Israel 
1980). 
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Table 4. Emission Optical Depth Maps Data 





Peak 


_max 
r 11.7 


Offset (Aa,A<5) 


< Til. 7 > 


r 24.5 


Offset (Aa,A<5) 


< T24.5 > 


Source 


Label 




(arcsec) 






(arcsec) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


11.11198-0.39795 


a 
b 


0.070 


(+1.4,-0.6) 


0.0055 


0.70 
0.60 


(+1.6,-0.8) 
(-0.2, +0.4) 


0.21 


11.94545-0.03634 


a 


0.0065 


(+0.8,+2.6) 


0.0013 










b 


0.0022 


(-1.2.+1.0) 




0.33 


(-0.3, +1.3) 


0.24 




c 


0.0008 


(+2.4,-1.7) 












d 








0.30 


(+2.2, +1.2) 




18.71179+0.00085 




0.003 


(+0.2.+0.8) 


0.0012 








19.75611-0.12775 


a 
b 


0.029 


(-0.2,+0.2) 


0.0053 


0.60 
0.47 


(+0.1,-0.1) 
(-0.7,-1.3) 


0.18 


21.38654-0.25346 




0.011 


(+0.7.+0.8) 


0.00066 








25.39918-0.14081 


a 
b 








0.80 
0.60 


(-1.6, +2. 3) 
(-2. 4, +1.1) 


0.18 




c 


0.019 


(+1.0,-3.8) 


0.0044 


0.57 


(0.0,-3.3) 






d 


0.014 


(+0.5,0.0) 












e 








0.56 


(+1.9,-0.8) 






f 








0.31 


(+3.8,-1.7) 






g 








0.24 


(+3.3,-3.7) 






h 








0.19 


(+4.2,-3.1) 




25.80211-0.15640 


a 


0.039 


(-0.1, +0.2) 


0.012 


0.33 


(0.0,+0.2) 


0.087 




b 


0.037 


(-0.8, +0.7) 




0.35 


(-0.8, +0.5) 




27.18725-0.08095 




0.005 


(+0.3,-0.2) 


0.0015 


> 1 


(+0.6,0.0) 


> 1 


28.28875-0.36359 


a 


0.12 


(+0.2,-1.6) 


0.033 


0.52 


(-0.2,-1.9) 


0.17 




b 


0.10 


(-0.3.+3.2) 




0.57 


(-0.2, +2. 4) 






c 


0.013 


(+1.2,-4.7) 












d 








0.49 


(+1.0,-1.4) 






e 








0.30 


(-3.7, +0.8) 




30.04343-0.14200 




0.0015 


(-0.2,0.0) 


0.00026 








30.86744+0.11493 




0.0034 


(-0.4, +0.3) 


0.00080 








30.58991-0.04231 




0.00028 


(+0.3,-0.6) 


0.00017 








30.66808-0.33134 




0.0097 


(+0.3,-0.3) 


0.0018 








33.91585+0.11111 


a 
b 
c 
d 
e 


0.012 
0.011 
0.008 


(+0.1,+2.7) 
(-l.l,+0.4) 
(+2.5,+4.5) 


0.0030 


0.08 
0.07 


(+2.1,0.0) 
(+1.1,+1.9) 


0.054 


33.81104-0.18582 


a 
b 


0.014 
0.0026 


(0.0,-1.9) 
(+0.1,-0.5) 


0.0030 


0.53 


(-0.5,-1.8) 


0.24 


35.46832+0.13984 


a 
b 

c 


0.038 
0.0061 
0.0055 


(-2.3.+0.3) 
(+2.4,+0.1) 
(+2.1, +2.8) 


0.0044 








37.87411-0.39866 


a 
b 


0.0075 
0.0059 


(+0.2,-0.1) 
(-2.7.+1.7) 


0.0015 










c 


0.0042 


(+0.2,-1.9) 




> 1 


(0.0,-1.4) 


0.15 




d 


0.0020 


(-3.8, +0.1) 
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Table 4 — Continued 



Peak 


max 
T 11.7 


Offset (Aa,A<5) 


< Til. 7 > 


r 24.5 


Offset (Aa,A<5) 


< T24.5 > 


Source Label 




(arcsec) 






(arcsec) 




(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


e 


0.0020 


(-4.6.+1.4) 










f 








0.45 


(-2.6,0.0) 




111.28293-0.66355 a 


0.023 


(-0.8,-1.5) 


0.0041 








b 


0.013 


(+0.5.+2.1) 










c 


0.008 


(+0.8,0.0) 











Note. — Optical depth peaks, the positions of the peaks, and optical depth means for our sources. The listed values 
are from our best alignment of the flux density maps. The typical Icr uncertainties due to a possible misalignment are 
200% for the peak optical depth and 50% for the mean optical depth. The dominant contribution to the uncertainties 
is from uncertainties in aligning the flux density maps. Some of the fields contain multiple optical depth peaks, which 
are labeled 'a', 'b', 'c', etc., and are indicated in the maps in Figures 4-31. 
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Table 5 lists the visible extinctions associated with the emitting dust. These were derived 
from the relation Ay 1 - 7 = 26.88 t u . 7 , or when possible from Ay 1 - 5 = 57.82 r 2 4. 5 (Draine 2003), 
where we use the average optical depths in columns (5) and (8) of Table 4, respectively. 
Many of our sources have Ay- 7 that are low enough to imply that many of the optical 
photons leak out of the region occupied by hot dust. The Ay 1 ' 5 values we derive for some 
of our sources are high enough to absorb all optical photons, thus being a better measure 
of the total extinction. However, the MIR colors and the MIR-derived spectral types of our 
sources indicate that at least for some them there is an additional extinction component that 
must come from cooler dust local to the source or on the line-of-sight (see §3.5 and §4.2). If 
dust and gas are well mixed, it is possible to derive the column density associated with the 
emitting dust by N H ~ 1.87- 10 21 A V [cm" 2 ] (Draine 2003). The column densities iV 1 /- 7 and 
iV^ 4 - 5 are given in Table 5 as well. See §3.5 for comparison with the column density of the 
ionized gas. 
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Table 5. Extinction and Column Densities of Emitting Dust 





411.7 


424.5 
A v 


A^ 1 - 7 x 10 21 


Nff- 5 x 10 21 


Source 


(mag) 


(mag) 


(cm -2 ) 


(cm -2 ) 


(1) 


(2) 


(3) 


(4) 


' (5) 


11.11198-0.39795 


0.15 


12 


0.28 


22 


11.94545-0.03634 


0.03 


14 


0.06 


26 


18.71179+0.00085 


0.03 




0.06 




19.75611-0.12775 


0.14 


10 


0.26 


19 




n D9 








25.39918-0.14081 


0.12 


10 


0.22 


19 


25.80211-0.15640 


0.32 


5 


0.60 


9 


27.18725-0.08095 


0.04 


> 58 


0.07 


> 108 


28.28875-0.36359 


0.89 


10 


1.66 


19 


30.04343-0.14200 


0.007 




0.01 




30.86744+0.11493 


0.02 




0.04 




30.58991-0.04231 


0.005 




0.009 




30.66808-0.33134 


0.05 




0.09 




33.91585+0.11111 


0.08 


3 


0.15 


6 


33.81104-0.18582 


0.08 


14 


0.15 


26 


35.46832+0.13984 


0.12 




0.22 




37.87411-0.39866 


0.04 


9 


0.07 


17 


111.28293-0.66355 


0.11 




0.21 





Note. — Visible extinctions and the corresponding column densities 
associated with the emitting dust for our sources. These figures were 
derived from the mean m.7 and T24.5. The typical la uncertainties are 
50%. 
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3.4. Morphologies and Distances 

Determining morphologies can be subjective and ambiguous, and with our sample we 
have only small numbers statistics, which makes it difficult to assess the statistical sig- 
nificance of the results. We use the morphological types that were defined by Wood & 
Churchwell (1989b) and revisited by Churchwell (2002), to determine the source morphology 
by eye. The relatively few morphological types of UC H II regions include: spherical or un- 
resolved, cometary, core-halo, shell, irregular or multiple-peaked, and bipolar. Table 6 lists 
our 18 sources with their morphology classes in column (2). We added comments explaining 
our choice of the morphological classes in column (3). 

Column (4) of Table 6 lists galactocentric distances (R ga i) which we derive using lit- 
erature line-of-sight velocities and the Galactic rotation curve of Rohlfs & Kreitschmann 
(1987). These velocities were determined using different methods and are based on various 
emission features (e.g., radio recombination lines; sub- mm CS molecular lines; OH, methanol 
and ammonia maser emission). Each entry in Table 6 has a reference to the source of the 
velocity. In some cases, there is also a reference to the work that resolved the distance 
ambiguity. In the two cases where no direct velocity measurements are available in the 
literature (18.71179+0.00085 and 30.66808-0.33134), we used molecular emission-line data 
from the Galactic Ring Survey (GRS; Simon et al. 2001) 3 . The GRS is a survey of the 13 CO 
J = 1 — > transition at 110.201 GHz, which is known to trace embedded UC H II regions 
(Koplak et al. 2003). The measurements are taken with a relatively fine sampling of 22" 
along the coordinate range / = 18°-54° and \b\ < 1, with a velocity resolution of 0.25 km s -1 . 
Because of the smaller line widths of 13 CO compared to 12 CO, it is possible to avoid velocity 
crowding and establish accurate kinematic distances to the clouds. For sources within the 
GRS coverage area, the GRS velocities and the velocities from the literature are consistent 
within the uncertainties implied by the line widths. 

The rotation curve of Rohlfs & Kreitschmann (1987) spans a range of 0.1 to 19.4 kpc 
in galactocentric radius. We have transformed this rotation curve (with R@ = 7.9 kpc and 
6 Q = 184 km s _1 ) to a rotation curve with the more updated IAU values, R Q = 8.5 kpc and 
# = 220 km s -1 . We assume that all sources lie in the Galactic plane (b = 0°), as all of 
them have very small latitudes < 0.7°). The kinematic distance (Rkin) is calculated by 
inverting the formula 

R 2 gal = R%- 2R km R Q cos / + R\ in . (2) 



3 This publication makes use of molecular line data from the Boston University Galactic Ring Survey 
(GRS). The GRS is a joint project of Boston University and Five College Radio Astronomy Observatory, 
funded by the National Science Foundation under grants AST-9800334, AST-0098562, & AST-0100793. 
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Table 6. Distances and Morphological Data 





Source 
(1) 




Morph. 
Class 
(2) 


Comment 
(3) 


R eii. 7 „ m 

(Kpc) (arcsccj 
(4) (5) 


Dll.7)ii 

(pc) 
(6) 




©5 GHz 
(arcsccj 
(7) 


D & GH: 
(PC) 

(8) 




11 


11198-0 


39795 


bipolar 


Extensions to N, SE 


17 ± l 1 


1.5 


0.12 ± 


01 


2.4 


0.20 ± 


.01 


11 


.94545-0 


03634 


shell 


Core at SE; Arc at NW 


4.2 ± 0.1 2 ' 3 


7.0 a 


0.143 ± 


004 a 


4.3 


0.088 ± 


.003 


18 


.71179+0. 


00085 


spherical 




4±1 4 


1.1 


0.021 ± 


005 


1.1 


0.021 ± 


.006 












12 ± l 4 




0.064 ± 


007 




0.064 ± 


.008 


19 


.75611-0. 


12775 


spherical 




4.2 ± 0.3 1 


0.8 


0.016 ± 


002 


1.8 


0.037 ± 


.003 












11.9 ± 0.3 1 




0.046 ± 


005 




0.104 ± 


006 


21 


.38654-0. 


25346 


comctary 


OfT-ccntcrcd core at SW 


5.6 ± 0.4 5 


1.3 


0.035 ± 


003 


2.3 


0.062 ± 


.005 












10.3 ± 0.3 5 




0.065 ± 


005 




0.115 ± 


.006 


25 


.39918-0 


14081 


bipolar 


Extensions to SSE and NNW 


9.8 ± 0.3 6 


4.0 a 


0.190 ± 


007 a 


3.7 


0.176 ± 


007 


25 


.80211-0. 


15640 


core-halo 


Extended halo 


5.5 ± 0.2 5 


1.2 


0.032 ± 


003 


1.3 


0.035 ± 


003 












9.8 ± 0.2 5 




0.057 ± 


004 




0.062 ± 


.005 


27 


18725-0 


08095 


irregular 


Weaker secondary peak at NE 


1.9 ± 0.3 1 


0.7 


0.006 ± 


001 


2.3 


0.021 ± 


.003 












13.3 ± 0.3 1 




0.045 ± 


006 




0.148 ± 


.007 


28 


.28875-0 


36359 


shell 


Shell with 2 peaks (17.65/iim) 


3.3 ± 0.1 1 ' 7 


5.0 a 


0.080 ± 


003 a 


4.1 


0.066 ± 


.003 


30 


.04343-0 


14200 


spherical 




6.1 ± 0.2 8 


0.6 


0.018 ± 


003 


3.2 


0.095 ± 


.004 












8.6 ± 0.2 8 




0.025 ± 


004 




0.133 ± 


.005 


30 


.86744+0. 


11493 


spherical 




12.2 ± 0.4 1 ' 9 


0.8 


0.047 ± 


005 


2.4 


0.142 ± 


.008 


30 


.58991-0 


04231 


spherical 




2.5 ± 0.8 10 ' 11 


0.8 


0.010 ± 


003 


2.5 


0.03 ± 


.01 


30 


.66808-0. 


33134 


spherical 




5.6 ± 0.3 4 


0.8 


0.022 ± 


003 


1.5 


0.041 ± 


004 












9.1 ± 0.3 4 




0.035 ± 


004 




0.066 ± 


.005 


33 


.91585+0. 


11111 


bipolar 


Extensions to NE, SW 


8.1 ± 0.1 12 ' 13 


3.0 a 


0.118 ± 


004 a 


3.3 


0.129 ± 


.004 


33 


.81104-0. 


18582 


irregular 


2 peaks aligned N-S 


11.2 ± 0.3 1,14 


3.0 a 


0.163 ± 


006 a 


1.1 


0.060 ± 


.006 


35 


.46832+0. 


13984 


bipolar 


Extension to E, W 


5.1 ± 0.3 1 


4.0 a 


0.100 ± 


006 a 


5.0 


0.124 ± 


.008 












8.8 ± 0.4 1 




0.171 ± 


009 a 




0.21 ± 


01 


37 


.87411-0. 


39866 


irregular 


Sub-peaks at S,NW 


9.5 ± 0.6 1 ' 15 


5.0 a 


0.23 ± 


02 a 


2.0 


0.092 ± 


.007 


111.28293-0.66355 


irregular 


Multi-peaks; filaments 


4.3 ± 0.5 1 ' 6 


2.0 a 


0.042 ± 


005 a 


4.7 


0.10 ± 


.01 



Note. — Morphologies (after Churchwcll 2002), comments on the choice of morphological class, kinematic distances, and 11.7 fim and 5 GHz 
diameters (Givcon ct al. 2005a). Multiple lines per object correspond to the near and far kinematic distances, in cases the distance ambiguity 
could not be resolved. 

a An upper limit on the FWHM since this object has an extended irregular or multi-peaked morphology. 

Note. — References: 1 - Bronfman, Nyman, &; May (1996). 2 - Braz 8z Sivagnanam (1987). 3 - Simpson Sz Rubin (1990). 4 - GRS (sec 
text). 5 - Szymczak, Hrynek, 8z Kus (2000). 6 - Churchwell, Walmsley, Sz Cesaroni (1990). 7 - Kurtz, Churchwcll, Sz Wood (1994). 8 - Lockman 
(1989). 9 - Palagi et al. (1993). 10 - Szymczak Sz Gerard (2004). 11 - Woutcrloot, Brand, Sz Ficglc (1993). 12 - Wink, Wilson, Sz Bicging (1983) 
13 - Wood Sz Churchwcll (1989a). 14 - Watson ct al. (2003). 15 - Afflerbach, Churchwcll, Sz Werner (1997). 
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For sources outside the solar circle, this gives a single solution, whereas inside the solar circle 
there are two solutions. The choice of distance in this case is based on different methods 
discussed in the references. 

MIR angular and physical sizes from this work are listed in columns (5)-(6). We resolve 
all sources. In most cases, we use the full width at half maximum (FWHM) as a measure 
for the source size, but in cases of extended irregular or multi-peaked sources, we estimate 
an upper limit on the FWHM. These cases are indicated by a footnote flag. Beam size (0.3" 
at 11.7 /im) was removed in quadrature to obtain the final size values. We also list radio 
angular (Giveon et al. 2005a) and physical sizes in columns (7)-(8) for comparison. 

3.5. Luminosities and Spectral Types 

MIR luminosities were calculated by integrating the Planck function from 1 to 1000 
(Urn using the dust color temperatures T17.65/11.7 and T24.5/17.65 and the corresponding optical 
depths. We employ the emissivity function (1 — e~ Tl/ ), with t v taken from Draine (2003), 
and assume isotropic emission into 4tt sr. Our MIR luminosities can be considered as lower 
limits to the bolometric luminosity, if we assume that most of the shorter wavelength flux 
from the central ionizing source is reprocessed by the dust. The results of §3.3 show that 
Ay derived from the 24.5 fxm data is sufficient to absorb all optical photons. That is usually 
not so for Ay derived from the 11.7 /im data. An anisotropic dust distribution, very thick 
dust that reprocesses radiation into the FIR regime (either locally or on the line-of-sight), 
and the reprocessing of radiation by the ionized gas in these sources are the main reasons 
for underestimating the bolometric luminosity, even at 24.5 /im. 

For each source, the Draine (2003) extinction curve was normalized to have the measured 
Tn. 7 or T24.5 of the source, when calculating £17.65/11.7 or £24.5/17.65; respectively. For this cal- 
culation we derive the total emission optical depths, T\\ .7 and r 2 4.5, for each source, based on 
the total flux densities at the corresponding bands. All sources except two (19.75611—0.12775 
and 27.18725—0.08095), have r v < 1 at 11.7 and 24.5 /im. These two sources are optically 
thick at 24.5 /im, so we estimate their luminosities by a blackbody with the corresponding 
temperature. 

We derive spectral types based on MIR flux densities and dynamical non-LTE model 
atmospheres of hot stars from Sternberg, Hoffmann, & Pauldrach (2003), assuming these 
bolometric luminosities come from a single star and that the MIR luminosity is a good 
estimator of the bolometric luminosity. For cooler stars (later than B0. 5) we use the tables 
of Doyon (1990), which are based on the Kurucz static LTE atmospheres (Kurucz 1992). MIR 
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luminosities and spectral types are listed in Table 7, columns (2)-(3) and (9), respectively. 
If the bolometric luminosity heating the dust comes from a cluster of stars, the MIR-derived 
spectral type of the ionizing star, which is typically the most massive star in the cluster, will 
be later than the listed type. 

An independent way to estimate stellar types is possible using our radio data, assuming 
they are dominated by free-free emission. Radio emission does not suffer from extinction at 
all, and thus allows us to test our assumptions regarding the MIR emission. We calculate 
the radio properties similarly to Wood & Churchwell (1989a) and list them in Table 7. For 
ionized gas, the brightness temperature at the Rayleigh- Jeans approximation at a frequency 
v is 

r - = i m L « 3 » 

where S v is the integrated flux density in milliJanskys and Q is the radio source solid angle 
in sr. 

The free-free optical depth, listed in column (4) of Table 7, is then calculated from 

T bu = (1 - e"-) T e , (4) 

assuming an electron temperature T e = 10 4 K with a 50% uncertainty. One exception for 
that is the source 37.87411—0.39866, which has a higher brightness temperature. In that 
case we assumed T e = 1.5 • 10 4 K. The emission measure, listed in column (5), is given by 

EM = — [pc cm- 6 ], (5) 

8.235 x 10- 2 a v T e -1 - 35 v G H y ' 

where a v is a unity order constant that equals 0.9938 for T e = 10 4 K and v — 5 GHz 
(Mezger & Henderson 1967). We can estimate the electron density by the rms value (Wood 
& Churchwell 1989a), 



10 6 pc cm" 6 / \D 

where D is the source diameter in pc. The densities are listed in column (6), and the column 
densities, N l ^ n derived from the densities are listed in column (7). Assuming a homogeneous, 
ionization-bounded H II region, in which none of the hydrogen ionizing photons ihv > 13.6 
eV) is absorbed by dust, we use the relation 

A"L y c = u B j n 2 e dV = a B EM J dtt R\ in (7) 



Table 7. Derived Physical Properties 





1 °S L 1T. 65/11. 7 a 


lo g L 24.5/17.65 a 






EMxl0 b 


n e X 10* 


N ton x 


log N Lyc 


MIR Spectral 


Radio Spectral 


Source 


(£©) 




(£©) 


T"5 GHz 




(pc cm b ) 


(cm" 3 ) 


(cm" 2 ) 




Typc b 


Type 


(i) 


(2) 




(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


11.11198 — 0.39795 


4 3 ± 


2 


5.2 ± 0.3 


0.15 ± 


05 


14 ± 5 


0.8 ± O.l 


2.5 ± 0.3 


48.5 ± 0.1 


B1-B0.5 / B0— 06.5 


08.5-09 


11.94545 — 0.03634 


3.4 ± 


1 


4.3 ± 0.2 


0.19 ± 


07 


17 ± 6 


1.4 ± 0.2 


1.9 ± 0.3 


47.9 ± 0.1 


B3— B2.5 / Bl 


BO. 5— BO 


18.71179+0.00085 








0.9 ± 


5 


80 ± 50 


6 ± 2 


1.9 ± 0.9 


47.3 ± 0.3 


B8— B5 


BO. 5 




3 5 + 












4 ± l 


4 ± l 


48.3 ± 0.2 


B4-B1.5 


09.5-09 


19.75611 — 0.12775 






4.0 ± 0.3 C 


0.07 ± 


02 


6 ± 2 


1.3 ± 0.2 


0.7 ± 0.1 


46.7 ± 0.1 


B4-B2.5 / B2-B1 


BO. 5 




4 2 + 




4.9 ± 0.3 C 








0.8 ± 0.1 


1.3 ± 0.2 


47.6 ± 0.1 


B1.5— B0. 5 / B0. 5— 08.5 


BO. 5 


21.38654 — 0.25346 


3 2 + 


2 




0.15 ± 


05 


14 ± 5 


1.5 ± 0.3 


1.4 ± 0.3 


47.5 ± 0.1 


B5-B2.5 


B0.5-B0 




3.7 ± 


2 










I . I ± 0.2 


2.0 ± 0.4 


48.0 ±0.1 


B2 


BO 


25.39918—0.14081 


4.3 ± 


1 


5.2 ± 0.2 


0.5 ± 


2 


50 ± 20 


1. 7 ± 0.3 


4.6 ± 0.8 


48.9 ±0.1 


B1-B0.5 / 09-07 


08— 07 


25.80211 — 0. 15640 


3.5 ± 


2 


4.2 ± 0.3 


0. 16 ± 


06 


14 ± 5 


2.0 ± 0.4 


I . I ± 0.2 


47.0 ±0.1 


B3— rS2 / B1.5— rSO.5 


BO. 5 




4.0 ± 


2 


4.7 ± 0.3 








1.5 ± 0.3 


1.4 ± 0.3 


47.5 ± 0.1 


B2-B1 / B0.5-O9.5 


BO. 5 


27.18725-0.08095 


1.9 ± 


4 


2.1 ± 0.2 C 


0.019 ± 


007 


1.7 ± 0.6 


0.9 ± 0.2 


0.29 ± 0.08 


45.6 ± 0.2 


A1-B8 / B9.5-B8 


B1.5-B1 




3.6 ± 


4 


3.8 ± 0.2 C 








0.34 ± 0.06 


0.8 ± 0.1 


47.3 ± 0.1 


B3-B1.5 / B2.5-B1.5 


BO. 5 


28.28875-0.36359 


4.2 ± 


1 


4.7 ± 0.2 


0.25 ± 


09 


23 ± 8 


1.8 ± 0.3 


1. 8 ± 0.3 


47.7 ± 0.1 


Bl / B0.5-B0 


BO. 5 


30.04343-0.14200 


3.7 ± 


7 




0.0014 ± 


0005 


0.13 ± 0.05 


0.12 ± 0.02 


0.18 ± 0.03 


45.8 ± 0.1 


B5-B0.5 


B1.5-B1 




4.0 ± 


7 










0.10 ± 0.02 


0.21 ± 0.04 


46.1 ± 0.1 


B3-B0.5 


Bl 


30.86744+0.11493 


4.1 ± 


4 




0.4 ± 


1 


36 ± 9 


1.6 ± 0.2 


3.5 ± 0.5 


48.6 ± 0.1 


B1.5-B0.5 


09-08 


30.58991-0.04231 


1.7 ± 


5 




0.06 ± 


02 


5 ± 2 


1.3 ± 0.3 


0.6 ± 0.2 


46.4 ± 0.2 


A4-B8 


B1-B0.5 


30.66808-0.33134 


3.2 ± 


3 




0.7 ± 


3 


60 ± 30 


4 ± l 


2.5 ± 0.7 


47.8 ± 0.2 


B5-B2.5 


B0.5-B0 




3.6 ± 


3 










3.0 ± 0.8 


3.1 ± 0.8 


48.2 ± 0.2 


B3-B1.5 


B0-O9 


33.91585+0.11111 


3.9 ± 


1 


4.7 ± 0.3 


0.3 ± 


1 


27 ± 9 


1.4 ± 0.2 


2.8 ± 0.4 


48.4 ± 0.1 


B2-B1.5 / B0.5-O9.5 


09.5-09 


33.81104-0.18582 


4.2 ± 


1 


4.8 ± 0.3 


0.6 ± 


3 


50 ± 30 


2.9 ± 0.9 


2.7 ± 0.9 


48.0 ± 0.2 


Bl / B0.5-O9 


B0.5-O9.5 


35.46832+0.13984 


4.0 ± 


1 




0.10 ± 


03 


9 ± 3 


0.8 ± 0.1 


1.5 ± 0.2 


47.9 ± 0.1 


B1.5 


B0.5-B0 




4.5 ± 


1 










0.6 ± 0.1 


1.9 ± 0.3 


48.4 ± 0.1 


B0. 5 


09.5-09 


37.87411-0.39866 


4.5 ± 


1 


5.1 ± 0.3 


1.3 ± 


7 d 


200 ± 100 


5 ± I 


7 ± 2 


49.0 ± 0.2 


BO. 5 / B0.5-O7 


08-06.5 


111.28293-0.66355 


3.6 ± 


2 




0.6 ± 


3 e 


4 ± 2 


0.6 ± 0.2 


0.9 ± 0.3 


47.3 ± 0.2 


B2.5-B2 


BO. 5 



Cn 
CO 



Note. — MIR luminosities, free-free optical depths at 5 GHz (or at 1.4 GHz for sources with no 5 GHz detection), emission measures, electron densities, column densities, Lyman continuum 
photon rates, and MIR and radio-derived spectral types for our sample. The MIR. luminosities and spectral types should be considered to be lower limits (see text). Multiple lines per object 
correspond to the near and far kinematic distances, as in Table 6. 

d Whcn a source was observed with all three bands, luminosities are derived for the two color temperatures calculated. 

b Spcctral types arc derived for the two luminosities calculated - L 17 7 / L 2 4 5/17 65- 

c Optically thick cases where luminosities were estimated using a blackbody spectrum. 

d T e = 15 000 K was assumed instead of 10 000. 

e This is actually Ti 4 GHz- 
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to calculate the Lyman continuum photon rate A^yc, 

^ = L7 x 10 « (* m \ (_°_a (**y ( 8, 

y \10 6 pc cm~ V Varcsec 2 J \kpc J 1 J w 

which is listed in column (8). No assumption was made here regarding the optical depth 
of the cloud. Assuming that these ionizing photons come from a single star, it is possible 
to determine the spectral type, again, using the stellar model atmospheres of Sternberg, 
Hoffmann, & Pauldrach (2003) for hot stars and of Doyon (1990) for cooler stars. The radio- 
derived spectral types are listed in column (10) of Table 7. It is unlikely that the central 
sources are star clusters, since our single-star estimates are already giving relatively late 
type stars. A central cluster would imply an even later type for the most massive star in the 
cluster. 

Tables 6 and 7 confirm that our sources, selected according to the combination of their 
radio and MIR continuum emissions are indeed UC H II regions, as they generally have very 
small physical sizes (D <, 0.1 pc), high densities (n e ^ 10 4 cm" 3 ), and have bright (EM ^> 10 6 
pc cm -6 , but in most cases > 10 7 pc cm -6 ) photoionized gas (Wood & Churchwell 1989a; 
Kurtz, Churchwell, & Wood 1994). 

Comparing dust-emitting column densities, A" 1 / 7 and A"^ 4 5 , derived from the LWS 
observations in this work (Table 5), and the ionized gas column density, N l ^ n , derived from 
our previous radio observations (White, Becker, & Helfand 2005; Giveon et al. 2005a) shows 
that Nff- 7 < Nfj n < A^ 4 5 . This supports the common assumption that the hot emitting 
dust is a relatively thin layer surrounding the Stromgren sphere of ionized gas, while the 
cooler dust that absorbs all of the optical photons consists of a much thicker column, further 
away from the central source. 

An interesting result that emerges from these calculations is that even according to the 
radio-derived ionizing luminosities, 50% of our objects are excited by B-stars and not by 
O-stars, and all our objects have derived types that are later than an 06 star. A possible 
explanation for missing the earliest O-stars might be a bias imposed by selecting the most 
compact radio sources (5 GHz diameters < 5"; see also §4.1). In 50% of the sources, the 
MIR-derived spectral types agree with the radio-derived spectral types. The 24.5 /zm-derived 
spectral types always show a better agreement with the radio ones compared to the 11.7 
//m-derived types. Two other sources (11.94545—0.03634 and 27.18725—0.08095) show clear 
deficiency in MIR photons (even at 24.5 /xm) to account for the radio-derived ionizing source. 
Both cases may be explained simply by the 24.5 /im luminosities underestimating the real 
luminosities due to the reasons given above. Alternatively, this may be explained by an 
Ay = 30 foreground dust exctinction, or by a comparable extinction by dust cooler than ~ 60 
K that is still local to the source. Seven sources that show a MIR-radio inconsistency in the 
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derived spectral types were not observed at 24.5 /xm. Their MIR luminosities, derived only 
using the 11.7 fim data, significantly underestimate the real luminosities. This inconsistency 
might be resolved by obtaining 24.5 /iin observations of these sources. 

On average, assuming the far distances for all sources that still have ambiguous dis- 
tances, the spectral type derived from £17.65/11.7 are later than the ones derived from £24.5/17.65 
- B2 (±2) compared to BO (±1, i.e., in the range Bl-09). Assuming the near kinematic 
distance for the ambiguous sources pushes the radio-derived spectral types to even later 
types, which are already quite late compared to what one expects from ionized compact HII 
regions. Non of the results changed significantly due to the choice of distance. The longer 
wavelength MIR-derived spectral types agree with the average radio-derived spectral type, 
which is also BO (±1). The later types derived from the 11.7 /iin luminosities (as late as 
A stars for 25.80211-0.15640 and 30.58991-0.04231), suggest that the derivation using the 
longer wavelengths may better match the real bolometric luminosities. The consistency of 
the spectral types derived from £24.5/17.65 and from the radio suggest that the contribution 
from plausible lower-luminosity non- ionizing companion stars in a cluster is not significant. 
This better consistency seems to be due to the 24.5 /xm emission coming from a thicker layer 
of dust compared to the 11.7 /itn - better representing the entire surrounding cloud - as 
is evident from our calculated emission optical depths (§3.3). Another reason for underes- 
timating the luminosities with £17.65/11.7 may be absorption by the broad 9.7 /im silicate 
feature going into the 11.7 /im band, and by the 18 /im silicate feature. 



4. Discussion 
4.1. Correlations 

We have tested for the existence of correlations among all measured quantities described 
above, both averages and peaks. We chose Spearman's rank-correlation coefficient over 
Pearson's correlation coefficient because it tests for a general monotonic relation rather than 
only a linear one. However, in cases of bimodality Spearman's coefficient may fail to detect 
correlations. 

The strongest correlation that emerges from our data is between the MIR and the radio 
integrated flux densities. Figure 32 shows the 24.5 /im flux densities as a function of the 5 
GHz flux densities with our best fit. 
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The correlation coefficient is r s = 0.94 with probability P r = 5 ■ 10~ 5 of being random. 
This strong correlation reflects the tight connection between the ionizing luminosity of the 
central star or stars, and the star's total luminosity reprocessed by the surrounding dust 
into the IR. However, the correlations of the 11.7 and 17.65 /im bands are much weaker - 
r,s = 0.54 (P r = 0.02) and r$ = 0.55 (P r = 0.02), respectively. These correlations remain 
weaker than that of the 24.5 /im even when only the same 10 sources are considered - 
r s = 0.87 (P r = 0.001) and r s = 0.74 (P r = 0.014) for 11.7 and 17.65 /im, respectively. 
The fact that of the three bands, the 24.5 fxm shows the strongest correlation with the radio 
further supports our result from §3.5 that the 24.5 /im flux density is a better estimator of 
the total flux density reprocessed by the dust. The source 28.28875—0.36359 is exceptionally 
bright in the MIR - perhaps part of the MIR emission in this extended object is not related 
to the ionizing source. Fitting a linear relation with la error weighting, we obtain, 

i^mfJy] = (0.38 ± 0.05) • F 5GHz [mJy] + (45 ± 8). (9) 

The source 28.28875—0.36359 does not affect the fit significantly due to the large uncer- 
tainty in its MIR flux density. Increasing the errors to 3a increases the errors of the fit - 
-^24.5^m[Jy] = (0.4 ± 0.2) • F 5 GHz[mJy] + (50 ± 20) - but the slope remains significantly larger 
than zero. 

Another correlation we find is between the 11.7 /im size D11.7 as defined in §3.4 and 
the ionizing photon rate, Ni yc : r s = 0.73 (P r = 6 ■ 10~ 4 ). D11.7 and N[ yc have the fitted 
power-law relation 

\ogN lyc [S- 1 ] = (1.98 ± 0.09) log A1.7 M + (50.1 ± 0.1). (10) 

The physical MIR size show weaker correlations with the MIR luminosities: r s = 0.63 
(P r = 0.005) with L17.65/11.7 and r s = 0.61 (P r = 0.008) with L 2 4.5/i7.65- Even though the 24.5 
/im flux densities show a stronger correlation with the radio flux densities, and the spectral 
types derived from the 24.5 /im flux densities are more consistent with the radio-derived types 
compared to the types derived from the 11.7 /im flux densities, there is a smaller number of 
sources at 24.5 /im and the range of their sizes is more limited. This makes the probability of 
their correlation with Ni yc smaller than the ll^um-iV^ correlation (r s = 0.87, P r = 0.001), 
and the uncertainty in the fit larger: log A^ c [s _1 ] = (2.4 ± 0.2) \ogD u . 7 [pc] + (50.3 ± 0.2). 
Also, most of the sources observed at 24.5 /im are still upper limits (6 out of 10) compared to 
a smaller fraction at 11.7 /mi (8 out of 18). We plot Ni yc vs. D n 7 with their la uncertainties 
in Figure 33, together with the best fit and its la uncertainty. 

This relation in eq. 10 implies that dust at the same temperature will be further away 
from the central source as Ni yc increases. This is because as the ionizing and heating power 
of the central star(s) increases, the dust destruction radius increases and the surviving dust 
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is heated at larger distances. The true power-law index is probably larger than 2, since a 
few of the sizes are upper limits (see §3.4). When the upper limits are not included, it is 
2.8 ± 0.2. Garay & Lizano (1999) found a correlation between the radio luminosity and 
the radio size for a collection of samples comprising a few hundreds UC and compact H II 
regions. Paladini, Davies, & DeZotti (2004) find a similar trend in a sample of 250 H II 
regions. The ionizing photon rate Ni yc is proportional to the radio luminosity (see eqs. 3-7), 
but other variables, such as the temperature and the angular size, introduce an additional 
scatter in Ni yc . 

Another important implication of this result is that the more compact sources are ac- 
tually due to later spectral types rather than young age. This was suggested using dynamic 
considerations by e.g., De Pree, Rodriguez & Goss (1995), Garcia-Segura & Franco (1996), 
Xie et al. (1996), and Garay & Lizano (1999), and was known as the age problem of UC 
H II regions (Kim & Koo 2001): their number is about an order of magnitude greater than 
expected from other indicators of massive star-formation rate based on their dynamical age 
(Garay & Lizano 1999). If most of the UC H II regions are in fact B-stars and not O-stars, 
as suggested by the present work, the age problem might be resolved. This is why we might 
be missing the earliest O-stars: selecting the most compact radio sources (5 GHz diameters 
< 5") may have biased our sample towards later types. 



4.2. MIR Color Criteria 

Wood & Churchwell (1989b) suggested IR color criteria for selecting UC H II regions, 
based on IRAS data. The ~ 150 factor in angular resolution between IRAS and LWS have 
an impact on these criteria, since in some cases, a significant fraction of the flux density 
comes from diffuse emission or sources unrelated to the compact object, that could not be 
separated with the IRAS beam. 

We compare the T25/12 temperature derived from the i^(25/xm)/i^(12/im) ratio of IRAS 
(Wood & Churchwell 1989b) to that we obtain from LWS (T24.5/11.7), assuming a modified 
blackbody with the Draine (2003) emissivity law as described in §3.2. These temperatures 
are shown in Figure 34. We preferred to compare the temperatures and not the raw flux 
density ratios since they are easier to interpret. The lower limit on the temperature for 
UC H II regions (163 K) is overplotted as a dashed line. This limit was derived from the 
limit on the F„(25pim) / F„(12pim) ratio (=3.7) from Wood & Churchwell (1989b), taking 
into account the bandpass response of IRAS and the assumed emissivity law. There is a 
significant decrease in the limiting temperature between IRAS and LWS to ~125 K. This 
decrease may be explained in part by the steep slopes of modified blackbody spectra at 100- 
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200 K, and the different band widths of IRAS and LWS, making the flux density in the IRAS 
12 fim band (8.0-14.5 fim band width) be dominated by emission at 14.5 fim - outside the 
LWS 11.7 yum band. However, simulating modified blackbody spectra at the temperature 
range 100-200 K show that for 2 out of the 10 sources with LWS T24.5/11.7 temperatures 
(11.94545—0.03634 and 27.18725—0.08095), the bandpass difference alone cannot explain 
the lower LWS temperatures. 

We conclude that at least in these two cases, this decrease is a result of the improved 
angular resolution: the larger beam of IRAS averaged out many lines-of-sights with both 
low and high extinctions, but the higher resolution of LWS means looking at the peak of the 
column density. Extinction local to the UC H II region may be caused, for example, by a 
broad silicate absorption feature centered at 9.7 /im, reaching into the 11.7 [im bandpass. A 
larger extinction at 11.7 /im compared to 24.5 /im will make the H II region appear cooler 
in the higher resolution observations. 

The original location of the limiting line in Figure 34 could have let in contamination 
from other source populations, which can explain the broader latitude distribution found 
using the lower resolution IRAS catalog (Wood & Churchwell 1989b) compared to works 
incorporating the better resolution of the MSX catalog (Giveon et al. 2005a). Such a change 
should be established using a larger sample of UC H II regions. A change in the Wood & 
Churchwell color criteria for selecting UC H II regions, will change their estimation of the 
formation rate of O stars in the Galaxy, and may lead to alleviation of the age problem of 
UC H II regions (§4.1). 

5. Conclusions 

We have presented high-resolution observations of a sub-sample of UC H II regions 
candidates drawn from the sample of Giveon et al. (2005a) in the first Galactic quadrant. We 
have presented sub-arcsecond resolution flux density maps, dust temperature and emission 
optical depth maps of these sources, and have studied the relations between their central 
stellar objects and the dust properties. The improved angular resolution of our observations 
(x50 compared to the MSX; xl50 compared to IRAS) lead to better pinpointing of the 
UC H II regions within the MIR sources, and thus, to better flux density and morphological 
determination. Some of the sources are still compact, but the peaks of all sources are resolved. 

Our main findings are: 

I. Half of our objects are excited by B-stars and not by O-stars, and all our objects have 
derived types that are later than an 06 star, even for radio-derived luminosities and 
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assuming far kinematic distances when ambiguous. Only two sources (10%) show a 
significant inconsistecy between the MIR and radio-derived spectral type, which may 
be explained by the MIR luminosities (I^s/is) strongly underestimating the bolomet- 
ric luminosities, or by foreground extinction of the MIR wavelengths (equivalent to 
Ay = 30). In general, approximating the bolometric luminosities using the longer 
wavelengths (17.65 and 24.5 /xm) is better compared to the shorter wavelengths (11.7 
and 17.65 fim). The average spectral type derived from L25/18 is BO (±1, i.e., in the 
range Bl-09) - the same as the radio-derived average spectral type. 

2. The majority of our sources are optically thin at MIR wavelengths, with a significantly 
optically thinner {r^f < 0.07 for all sources) hot dust (150-250 K), compared to an 
optically thicker {t™£ > 0.08 for all sources and > 1 in two cases) component of warm 
dust (60-110 K). Most of the optical photons are absorbed by the warm dust. 

3. The 24.5 /im flux densities of our sources are correlated with the 6cm radio flux densi- 
ties. This relation is fitted by F 24 .5^ m [Jy] = (0.38 ± 0.05) • F 5G Hz[mJy] + (45 ± 8). The 
correlations of the 11.7 and 17.65 fim bands are weaker, further supporting our result 
that the 24.5 /im flux density is a better estimator of the total flux density reprocessed 
by the dust. 

4. The MIR sizes are correlated with the source ionizing photon rates, confirming the 
picture of dust cocoons enveloping Stromgren spheres of ionized gas. We find the 
relation \ogN lyc [s" 1 ] = (1.98 ± 0.09) log D u , 7 [pc] - (50.1 ± 0.1). This correlation 
suggests that the more compact sources are actually due to later spectral types rather 
than young age, which may lead to alleviation of the age problem of UC H II regions. 
A possible explanation for missing the earliest O-stars might be a bias imposed by 
selecting the most compact radio sources (5 GHz diameters < 5") . 

5. The new flux density estimations lead on average to redder MIR colors or cooler temper- 
ature for our sources compared to IRAS observations. For at least two of our sources, 
the reddening effect is too big to be exclusively explained by the different bandpass 
response of IRAS and LWS. We conclude that in these cases the significantly smaller 
beam of LWS focuses on larger dust columns towards the more accurate positions of 
the UC H II regions, while the larger beam of IRAS averaged out many lines-of-sight. 
Our sources are thus shifted to redder colors in the Wood & Churchwell (1989b) color 
diagram that is used for selecting UC H II regions based on their MIR colors. This 
may account for some contaminations in selecting UC H II regions using the IRAS 
criteria, and may lead to a change in the estimated star-formation when larger samples 
are observed at sub-arcsecond resolution. 
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6. The MIR maps show a good overall match in shape when compared to existing radio 
images with comparable resolution. 

The significant improvement in angular resolution and the ability to pinpoint the UC 
H II regions lead us to conclude that such an analysis applied to a larger sample would 
improve our understanding of star-formation in the Milky Way: what are the effects of 
including a possible large population of B-stars, and the reddening effect we discovered on 
color-selection criteria of UC H II regions in both the Galactic and extragalactic objects? 
Our observations had a 25% rate of non-detections with LWS of sources with high MIR- 
radio matching probabilities based on the MSX and the VLA data. We suspect the reason is 
the large difference in angular resolution between LWS and MSX, combined with the small 
field of view of LWS. This either caused a false match between the IR source and a radio 
source, or even if the match is real, the IR source might be a low surface brightness source, 
below the detection limit of LWS. The GLIMPSE survey is most suitable for confirming 
the MIR-radio associations of sources in our original sample (Giveon et al. 2005a), since 
it has higher resolution than the MSX. However, it is limited in wavelength. Additional 
observations which will compensate for that deficiency are thus needed. 
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Fig. 32. — Flux density at 24.5 jim. as a function of the 5 GHz flux density with la errors. 
The horizontal error bars are smaller than the plotted circles. The linear fit given in eq. 9 
is overplotted as a dashed line. The solid lines mark the la uncertainty around this fit. 
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Fig. 33. — The ionizing photon rate as a function of the 11.7 /im physical size with la 
errors. Source MIR size increases with the ionizing rate of the central star. The fit to this 
relationship, expressed by eq. 10 is overplotted as a dashed line. The solid lines mark the 
la uncertainty around this fit. 
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Fig. 34. — A comparison between T 25 / 12 derived from the i^(25/xm)/i^(12/xm) ratios of 
sources in our sample (when available) as measured by IRAS (empty circles) and by LWS 
(filled circles), with la error bars. Arrows indicate limits on the ratios. The vertical dashed 
line is the lower limit of 163 K for UC H II regions according to Wood & Churchwell (1989b). 



